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Abstract 

It is highly desirable to integrate graphene into existing semiconductor technology, where the 
combined system is thermodynamically stable yet maintain a Dirac cone at the Fermi level. First- 
principles calculations reveal that a certain transition metal (TM) intercalated graphene/SiC(0001), 
such as the strongly-bound graphene/intercalated-Mn/SiC, could be such a system. Different from 
free-standing graphene, the hybridization between graphene and Mn/SiC leads to the formation 
of a dispersive Dirac cone of primarily TM d characters. The corresponding Dirac spectrum is 
still isotropic, and the transport behavior is nearly identical to that of free-standing graphene for 
a bias as large as 0.6 V, except that the Fermi velocity is half that of graphene. A simple model 
Hamiltonian is developed to qualitatively account for the physics of the transfer of the Dirac cone 
from a dispersive system (e.g., graphene) to an originally non-dispersive system (e.g., TM). 
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Graphene, a one-atom-thick carbon sheet, has sparked enormous interest during the past 
few years. Its peculiar band structure, with two linear dispersion bands crossing at the 
Fermi level (Ep) and other outstanding physical properties, make graphene an appealing 
system for both fundamental studies and modern technological applications UQ. Due to 
the two-dimensional nature of graphene, however, a substrate is required for most practical 
applications. The coupling with the substrate, however, leads to modification of the graphene 
morphology in the form of lattice corrugation, which could significantly affect its electronic 
properties MS . It is generally believed that to preserve the Dirac cone, graphene should 
only bind weakly to the substrate via, for example, van der Waals interactions. From a 
technological point of view, however, the stronger the binding, the better the thermal and 
mechanical stabilities of the combined system. It is, therefore, imperative to find Dirac cone 
structures that can also exist in strongly-bound systems with significantly larger binding 
energy. 

In searching for a solution, two issues must be addressed: (i) how to increase the binding 
without significantly increasing the corrugation of the graphene. For example, graphene/SiC 
(G/SiC) can be considered as a strongly-bound system due to the formation of two C- 
Si bonds out of eight carbon atoms per 2x2 graphene cell {9]. However, because only 
25% of the carbon atoms are involved in such a bonding, the net result is a heightened 
corrugation of the graphene. In this regard, a buffer layer with less directional bonding 
between graphene and SiC is highly desirable, (ii) A chemically active buffer layer may cause 
considerable charge transfer between graphene and substrate. Because graphene is strongly 
electronegative, most likely its Dirac cone will be filled with electrons. One therefore needs 
to engineer a new Dirac cone in the vicinity to replace that of graphene; in other words, 
to transfer the Dirac fermion of graphene to the surface of the substrate. For this to work, 
clearly, the surface of the substrate must possess at least doubly degenerate states: for 
example, p states at the symmetry points of the surface Brillouin zone, or transition metal 
(TM) d states, or rare-earth metal / states. 

Using first-principles calculations, we have systematically studied the 
graphene/intercalated TM/SiC(0001) systems (denoted as G/i-TM/SiC) with 3d TM 
elements. We find clearly-defined Dirac cone structures centered at the K point of the 
Brillouin zone as a consequence of the interaction among graphene, SiC s, p states and 
degenerate TM d states. In contrast to weakly-bound systems, Dirac states are primarily 
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located in the TM layer with diverse electronic properties. For example, G/i-Mn/SiC is a 
non-magnetic system with the Dirac point exactly at Ep\ G/i-Cr/SiC is a spin-polarized 
half Dirac fermion system in the sense that there is a Dirac point for majority-spin states 
0.2 eV above Ep but not for minority-spin states; G/i-Co/SiC and G/i-Ni/SiC show a 
gap at the Dirac point due to symmetry lowering, yet both maintain a good linear energy 
dispersion. A band coupling model is constructed to qualitatively explain the transfer of the 
Dirac fermion from graphene to TM in these coupled systems. The significant hybridization 
between graphene, SiC, and TM also results in noticeably stronger binding of graphene to 
the substrate with a binding energy Eb as high as 0.51 eV/C for G/i-Mn/SiC. Transport 
calculation reveals that the Dirac fermion of G/i-Mn/SiC behaves almost the same way as 
that of free-standing graphene, signifying the high potential of the G/i-TM/SiC systems in 
future electronic and spintronic applications. 

Our electronic structure calculations are performed within the local spin density approx- 
imation (LSDA) 10[ as implemented in the Vienna ab-initio simulation package (VASP) 
The electron-ion interaction is described by a projector augmented wave method 12 1 
with an energy cutoff of 400 eV. Integration over the Brillouin zone is carried out using the 
Monkhorst-Pack scheme with 12 x 12 x 1 Appoints. A x y/3R30° supercell is used for 
6/f-SiC, which can accommodate a 2 x 2 graphene cell, with one TM atom intercalated 
between the graphene and SiC substrate (see Fig. 1). The SiC substrate is modeled by a 
six SiC bilayer slab with hydrogen passivation of the bottom surface. In the simulation, the 
three bottommost SiC bilayers are fixed at their respective bulk positions whereas all other 
atomic positions are fully relaxed without any symmetry constraint until the residual forces 
are less than 0.01 eV/A. 



Quantum transport calculations are conducted using the ATK package [13J, which imple- 
ments nonequilibrium Green 's-f unction formalism within density-functional theory (DFT). 
The Eb of the graphene with the TM/SiC substrate (per carbon atom) is defined as: 
E h = (E G + E TM/SiC -£ G /i-TM/sic)/8, where E G , E TM/SiC and E G/i _ TM/SiC are the total en- 
ergies of isolated graphene, TM/SiC, and the combined G/i-TM/SiC systems, respectively. 
The positive value here means that the process is exothermic. 

In all the G/i-TM/SiC systems, it was found that the TMs bind to three surface Si 
atoms with a local C% v symmetry. In contrast, the positions of the TMs with respect to the 
hexagon of graphene vary. For example, Sc, Ti, V, Cr, Mn, and Fe reside directly below 
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FIG. 1: (Color online) Top view of the optimized geometries of the G/i-TM/SiC systems, (a) Sc, 
Ti, V, Cr, Mn and Fe, (b) Co, and (c) Ni. (d) Side view of the geometry corresponding to (a). 
Yellow, gray, and white balls are Si, C, and H, respectively, while others are TM elements. For 
clarity, graphene is shown as hexagonal mesh. Rhombus denotes the supercell. 

the center of the hexagon [see Fig. 1(a)], Co is about half way away from the center [see 
Fig. 1(b)], whereas Ni is on the edge directly below a C-C bond [see Fig. 1(c)]. These are 
distinctly different from the adsorption of TMs on graphene, where the low-energy site is 



always the center of the hexagon [14j . As compared to G/SiC, the corrugation of graphene 
is significantly reduced to practically negligible for Sc, Ti, V, Cr, Mn, and Fe [see Fig. 1(d)], 
and to only ~0.1 A for Co and Ni. The separation between the TM layer and graphene is 
~1.6 A in most cases except for Ni where the separation is 1.88 A. 

We characterize G/i-TM/SiC as strongly-bound graphene systems: local density approx- 



imation (LDA) calculation reveals that Ej> for graphene on TM/SiC > 0.24 eV/C 15], which 
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is at least three times the calculated adsorption energy of graphene on Pd(lll) using the 
same approach [5]. Graphene on Mn/SiC has the largest Eb = 0.51 eV/C. This value is 
even larger than that of graphene on SiC of 0.36 eV/C (9). It is known that LDA typically 
overestimates Eb, while generalized gradient approximation (GGA) typically underestimates 
Eb, although the relative energies will not change significantly. Our GGA 16) value for G/i- 
Mn/SiC is Eb = 0.35 eV/C, which can be contrasted to the GGA result for van der Waals 
binding — typically less than 0.1 eV/C. Therefore, the binding of graphene to TM/SiC is 
indeed reasonably strong. To experimentally realize the G/i-TM/SiC structures, one ap- 
proach is to make use of the fact that metal atoms can diffuse through the honeycomb of 



epitaxial graphene 17H19|. Another approach is to deposit a sub-monolayer TM on the SiC 



surfaces before covering them with graphene. Upon increasing the temperature, the TM 
atoms are expected to diffuse to the lowest-energy positions shown in Fig. 1. 

Figure 2 shows the calculated band structures for G/i-TM/SiC. The main feature of 
the bands can be characterized as the existence of two linear bands in the vicinity of Ep, 
highlighted by magenta lines. For Sc, Ti, V, Cr, Mn and Fe, the two bands cross at the 
K point, forming a Dirac cone. For Co and Ni, the Dirac cone splits, resulting in a small 
gap of approximately 0.15 eV. Going from Sc to Ni, the position of the Dirac point shifts 
from above to below Ep, indicating a conduction-type change from p to n. G/i-Mn/SiC is 
a special case where the Dirac point lies exactly at Ep. G/i-Co/SiC is also unique in that 
it has a direct gap of 0.26 eV at the M point. Although both TM-doped graphene 14j and 
TM-doped SiC [2Q] are magnetic systems, most of the systems we studied are non-magnetic. 
G/i-Cr/SiC is the only exception. A fully spin-polarized Dirac cone may exhibit negative 
refraction for one spin, while maintaining positive refraction for the other spin. When an 
electron beam travels through an interface between regions with different types of carriers, it 
can be deflected and, depending on the shape of the interface, focused. Combining these two 
effects may give rise to focused spin-polarized electron beams (spin lens) [2jJ. G/i-Cr/SiC 
would be such a system for experimental realization of spin lens. 

At first glance, the Dirac cone seems to originate from graphene. Two pieces of evidence, 
however, indicate that this cannot be the case. First, metals have significantly smaller 
workfunctions than graphene. From our calculations, the workfunctions of TM/SiC are all 
less than 3.3 eV, in contrast to ~4.6 eV for graphene {22]]. As a result, electrons are always 
transferred from TM/SiC to graphene. Charge density analysis reveals that the amount of 
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FIG. 2: (Color online) Band structures of the G/i-TM/SiC systems. Bands having a linear (or 
nearly linear) dispersion around the K point and near the Ep are highlighted by magenta except 
for spin-polarized G/i-Cr/SiC for which black solid lines denote majority-spin states and green 
dashed lines denote minority-spin states. Fermi levels are indicated by the red horizontal dashed 
lines at energy zero. 

charge transferred in G/i-TM/SiC is typically 0.5 e/TM. Thus, a Dirac cone that is primarily 



attributed to graphene cannot be above the Ep, as in the case of Sc, Ti, V and Cr. T 



qual.tat.vely dufeeut from decoupled graphene by iute r ca.atiu g H |22|, Ge M, Li Qll* 
or Au [19[ between graphene and SiC. Second, it has been shown that a strong coupling 



his is 



between graphene and TM may destroy the linear dispersion of graphene [5j. Therefore, 
there must be some yet unknown mechanism responsible for the observation of the Dirac 
cones in Fig. 2. 

As a prototypical example, let us examine G/i-Mn/SiC more closely. Previous work has 
shown that Mn binds weakly to graphene due to the high Mn 4s-to-3d electron promotion 
energy [ijj]. Therefore, it has to be the SiC substrate that alters the chemistry of Mn to 
increase its binding to graphene. Figure 3(a) shows schematically the interactions between 
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FIG. 3: (Color online) (a) Schematic drawing of the hybridization between Mn orbitals and the 
three dangling-bond orbitals of surface Si in Mn/SiC. Each arrow denotes one electron with up 
or down spin, while each black dot denotes two (spin-degenerate) electrons, (b) Local density of 
states (LDOS) of G/i-Mn/SiC. Zero is the Fermi level, (c) Schematic drawing of the hybridization 
between the Dirac cone of Mn/SiC and that of graphene: left is before and right is after. 

Mn and Si dangling-bond orbitals, derived from the calculated local density of states (LDOS) 
for Mn/SiC. Note that in the calculation, Mn is kept at the position before graphene removal. 
The interactions lead to six hybrid states, ax, e 2 (2), ei(2), and ai(s), respectively, with 
predominantly metal character: the ax singlet consists of Mn 4s and d z 2 character; the e 2 
doublet consists of Mn d xy and d x 2_ y 2 character; the ex doublet consists of Mn <i xz and d yz 
character; and the ai(s) singlet consists of Mn 4s character. Only the ax and e 2 states have 
noticeable contributions from Si dangling bond orbitals. Ten electrons, seven from Mn and 
three from Si dangling bonds, fully occupy the ax, e 2 and ex states, respectively, leaving the 
high- lying ai(s) state empty. Critically important is that both the <i-band e 2 and ex states 
become Dirac cones at the K point of the Brillouin zone, as schematically shown in Fig. 3(a). 
Therefore, the presence of SiC gives rise to two important effects: (1) forcing the donation 
of Mn 4s electrons to Mn 3d orbitals and (2) mediating the formation of (relatively-flat) 
<i-band Dirac cones. 
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Next, we consider the hybridization between graphene and Mn/SiC. Figure 3(b) plots 
the LDOS for G/i-Mn/SiC, which shows clearly the dominant Mn d-character of the Dirac 
cones in Fig. 2. As mentioned earlier, before hybridization both graphene and Mn/SiC 
have their Dirac cones. Because of the workfunction difference, the Dirac point of Mn/SiC 
[denoted as D^ in Fig. 3(c)] is above that of graphene [denoted as Dq]. Upon hybridization, 
the (originally nearly flat) Dirac cone of Mn/SiC is significantly broadened, whereas that of 
graphene is significantly narrowed, as deduced from the calculated band structure of G/i- 
Mn/SiC in Fig. 2. To qualitatively understand the effects of hybridization, we constructed 
a 4 x 4 matrix to mimic the interaction between two Dirac cones: 

/ -e + c G k a b \ 

— e — c G k b a 

(1) 

a* b* ctm^ 

\ b* a* -c TM k J 

where h is set to 1, A; is the wavevector measured from the K point, e > is the energy 
separation between the two Dirac points, cq and ctm are the speeds of Dirac fermions for 
graphene and TM, respectively, and a and b are coupling constants between states of the 
same and opposite directions of velocity. While the matrix can be numerically diagonalized, 
it is instructive to examine the solutions when a = and separately when 5 = 0. 

When a = and \b\ e, the hybridization results in a uniform reduction in the speeds 
of graphene and TM by a factor |6/e| 2 (cc + ctm)- When 6 = and \a\ <C e, on the other 
hand, the coupling results in a reduction in the speed of graphene by |a/e| 2 (cc — ctm) but 
an increase in the speed of TM by \a/e\ 2 (cQ — ctm)> as we have empirically observed. Our 
model fitting to the calculated results, without the small a assumption, shows that the 
dispersion for Mn in Fig. 2 corresponds roughly to b = and a = 3.5 eV for e = 1.5 eV, 
c G = 1-05 xlO 6 m/s, and c TM = 10 5 m/s. We note that the actual coupling may involve 
multi-bands and is hence more complex. Our fitting here is only to provide a physical 
picture. For a high crystal symmetry as the one in Fig. 1(a), it is conceivable that the 
coupling between waves of opposite traveling directions, namely, b, is small. However, for a 
lower crystal symmetry as in the case of Co in Fig. 1(b) or Ni in Fig. 1(c), however, this is 
highly unlikely. When neither a nor b is zero, our model suggests a gap at the Dirac point: 
Eg = ( v /e 2 + 4(a + 6) 2 - ^e 2 + 4(a - b) 2 )/2 w 4a6/Ve 2 + 4a 2 for b < a, e. The presence of 
a gap is consistent with the results of direct calculations for Co and Ni, showing a gap of 
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FIG. 4: (Color online) (a) Band structure around the Dirac point for G/i-Mn/SiC. Left: three- 
dimensional view and right: constant energy contours. Color represents vertical distance from 
zero energy surface, (b) The I-V curves for 11-AGNR (= armchair graphene nanoribbon with 11 
dimer rows) with and without the Mn/SiC support. Inset shows the atomic structure for transport 
calculations. 

about 0.15 eV. Graphene has about 8% lattice mismatch to SiC js]. Using a x Q\/3R30 o 
surface cell, which is nearly lattice matched to a 13 x 13 graphene cell, we find a gap of 
0.16 eV for G/i-Mn/SiC. The gap opening here can be attributed to an increased lattice 
corrugation similar to that of Co and Ni. Despite the gap opening, our fitting also reveals 
that the linear dispersion remains intact as long as the wavevector k (measured from the K 
point) is larger than 0.01 n/a. These small-gap G/i-TM/SiC systems may be used as planar 
linear-band semiconductors. 

Due to the symmetry of the d orbitals, one may expect the <i-state-derived Dirac cone 
to be highly anisotropic 25j. If so, group velocity may be significantly altered, at least in 
certain directions. The left panel of Fig. 4(a) shows a three-dimensional plot of the Dirac 
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cone for G/i-Mn/SiC. It reveals that the Dirac cone is as isotropic as that in graphene with 
a Fermi velocity of 0.5 x 10 6 m/s, which is almost half that in free-standing graphene. A 
linear dispersion is maintained for energies within ±0.3 eV from the Dirac point, which is 



more than half of that for graphene, ±0.5 eV [26j. The constant-energy contours of the 
valence band around the Dirac point is shown in the right panel of Fig. 4(a). It confirms 
that all the contours are concentric circles as are those in graphene. Figure 4(b) shows 
the calculated current- voltage (I-V) curves. It can be seen that the I-V characteristic of 
G/i-Mn/SiC under a moderate bias of less than ~0.6 V is also similar to that of graphene. 
Thus, a free-standing graphene-like transport can also be achieved on such strongly-bound 
systems. 

In summary, first-principles calculations reveal Dirac spectra in strongly-bound G/i- 
TM/SiC systems. The SiC substrate modifies the electronic structure of the TM d states to 
foster a strong coupling between TM and graphene. This also leads to a migration of the 
Dirac cone at the Fermi level from graphene to TM, which can be qualitatively understood 
via a simple model Hamiltonian. In the case of G/i-Mn/SiC, the Fermi level is at the Dirac 
point. Despite the d character of the states, the Dirac fermions behave the same way as those 
in graphene except for a half Fermi velocity. For other TM, our study also reveals interest- 
ing physical phenomena such as symmetry breaking, gap opening and spin-polarized Dirac 
cone near Ep. These results point to alternative and perhaps more realizable directions for 
graphene-based electronics and spintronics. 
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